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We analyse the dynamics of the Warsaw Stock Exchange index WIG at 
a daily time horizon before and after its well defined local maxima of the 
cusp-like shape decorated with oscillations. The rising and falling paths 
of the index peaks can be described by the Mittag-Leffler function super- 
posed with various types of oscillations. The latter is a solution of our 
model of index dynamics defined by the nonhomogeneous fractional relax- 
ation equation. This solution is a generalised analog of an exactly solvable 
model of viscoelastic materials. We found that the Warsaw Stock Exchange 
can be considered as an intermediate system lying between two complex 
ones, defined by short and long-time limits of the Mittag-Leffler function; 
these limits are given by the Kohlraush- Williams- Watts law for the initial 
times, and the power-law or the Nutting law for asymptotic time. Hence 
follows the corresponding short- and long-time power-law behaviour (differ- 
ent " universality classes" ) of the time-derivative of the logarithm of WIG 
which can in fact be viewed as the "finger print" of a dynamical critical 
phenomenon. 

PACS numbers: 05.45. Tp, 89.65.Gh, 89.75.-k 

1. Introduction 

It seems that there are many distinct analogies between the dynamics 
and/or stochastics of complex physical and economical or even social sys- 
tems UlllElilinilllZllHlEl- The methods and even algorithms that have 
been explored for description of physical phenomena become an effective 
background and inspiration for very productive methods used in analysis of 
economical data ^1 . 
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In this paper we study an emerging market and more precisely, the his- 
torical Warsaw Stock Exchange (WSE) index WIG at a daily time horizon 
at the closing; we think that its dynamics is typical for an emerging finan- 
cial market of small and moderate size. Our concept is to consider only well 
developed temporal (local) maxima of the cusp-like shape decorated with 
some oscillations (cf. peaks denoted by A, a, B, C in Fig^ which cover the 
greater part of the whole time series. Our goal is to describe the slowing 
down of rising and relaxation processes within these temporal maxima by 
assuming a retarded feedback as a principal effect (except the rising path 
within the first local peak which, in principle, is exponential; cf. Fig|2I). 
This feedback is a reminiscence of investors' activity stimulated mainly by 
their observations of the empirical data in the past. 



1.1. Inspiration 

Our analysis was inspired by the non-Debye or non-exponential, frac- 
tional relaxation processes observed, for example, in stress-strain relaxation 
present in viscoelastic materials |191 1131 114j . The most commonly used em- 
pirical decay function for handling non-Debye relaxation processes in com- 
plex systems are described either by a Kohlraush- Williams- Watts (KWW) 
or a stretched exponential decay function for short time-range t <C t: 

/(i)~exp(-(^)"), (1) 

where < a < 1, or by an asymptotic power-law (Nuttig law) for t ^ r: 

X 



. (2) 



time^ 



From expressions and © we obtain for the change of ln/(t) per unit 
din fit) r 1 fort«r 



dt It for t » T, 



(3) 



which defines the power-law limits of the dynamics, i.e. it defines (for t <C r) 
a universality class (characterized by dynamical, critical exponent 1 — a) 
which, in fact, can be viewed as the "finger print" of a dynamical critical 
phenomenon. For such a peak the derivative diverges at the extremal point 



^ The logarithm of the index or price is the quantity playing a fundamental role in 
financial analysis both of stochastic or deterministic types, e.g. its time-derivative is 
the instantaneous interest rate or return per unit time. 
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(i.e. at t = 0) which justifies the name "sharp peak" used often in the 
Uterature in this context (cf. [H] and refs. therein). For a sufficiently 
wide empirical window it would be possible to observe a transition from the 
KWW to the Nutting behaviour and from one kind of power-law to another, 
correspondingly. 

Note that the non-exponential relaxation introduces memory, i.e. the 
underlying fundamental processes are of non-Markovian type. It was shown 
that a natural way to incorporate memory effects is fractional calculus. 
The power-law kernel defining the fractional relaxation equation involves a 
particularly long memory. The function which plays a dominating role in 
fractional relaxation problems is indeed the Mittag-Leffler (ML) function 



which is a natural generalisation of the exponential one. This function in- 
terpolates between cases and and plays a central role in our analysis. 



Our phenomenological model of index dynamics consists of two stages: 

(i) Formulation of a linear ordinary differential equation of the first or- 
der with no feedback incorporated which describes evolution of an 
auxiliary, synthetic index only. 

(ii) A conjecture which transforms the above mentioned equation to a 
more general fractional form which already models the evolution of 
the empirical Warsaw Stock Exchange index. 

The transition from (i) to (ii) means that the system is changed from an 
unrealistic to a realistic, complex one where the retarded feedback plays an 
essential role. By using this model we are able to describe the well developed 
temporal maxima present in daily time series defined by WIG activity (cf. 
four maxima A, a, B and C shown in Fig^. 





2. The model 



2.1. Evolution of WIG 



The time-dependent value of WIG, X{t), can be decomposed into two 
components which are recorded and can be even electronically accessible for 
traders: 
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Fig. 1. The daily closing value of the Warsaw Stock Exchange index WIG (measured 
by conventional points (p)) from April 16, 1991 (the beginning of the Warsaw Stock 
Exchange activity) to the 30*'* of March 2006; the presented time series consists 
of empirical data points for 3382 trading days (td). The local maxima denoted as 
A, a, B and C are analysed further in the text. 



(a) The instantaneous offset U{t) =' Uj^{t) — U-{t) between the temporal 
(total) demand ?7+(t) > for stocks defining WIG and their (total) 
temporal supply U^{t) > 0. 

(b) The instantaneous volume trade V{t) min[U+{t), U-{t)] of stocks 
of the companies which consitute WIG. 

Hence, we can write the following instantaneous superposition 



where A{> 0) and B are referred to as coefficients of relocation. Note that 
paradoxally index X{t)'^ can change even when the volume of trade, V{t), 
vanishes, i.e. when the demand or supply vanishes. On the other hand, for 
vanishing U (i) (when the demand is balanced by the supply) volume trade, 
V{t), can be still nonvanishing which leads to the change of WIG'^. 

Note that X{t) is defined here relative to some reference level so it can be, in general, 
both positive and negative. 

Eq.Q defines an additive variant of our model though a multiplicative one is also 



X{t) = A-U{t) + B - Vit), 



(5) 



possible. 
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To consider the dynamics (evolution) of WIG, we complete eq.© by 
the differential one which exhibits an instantaneous, therefore unrealistic, 
feedback to the financial market, namely 

where coefficients C, D are rates while E is again a kind of relocation 
coefhcient. By combining eqs.© and © we eliminate the V{t) variable 
and obtain (after integration), 

X{t) - X(0) - sgn{C') rf ^ oA"'^(i) = -A' sgn{C) r^^ oDi'U{t) 

+ A'[Uit)-Um, (7) 

where the definition of an inverse derivative of the first order was used; the 
definition of its general n*'^ order version (for n = 1,2,3,...) has a useful 
form given by the Cauchy formula of repeated integration 

0A""/(*) = /* dtn-l ' dtn-2 . . . r f{to)dto 

Jo Jo Jo 

= -3- f dt'{t-tr-'f{t'). (8) 

i {n) Jo 

The combined coefficients 

which are valid for BE ^ 1; otherwise, instead of eq.Q we would obtain a 
very special one. 

The integral eq. ((7j) defines the model which is an analog of the Zener one 
for viscoelastic solids |17| in which the stress {U) - strain {X) relationship'* 
is given originally by the linear first order differential equation jl4j . Indeed, 
eq.(I7I) is the one which we generalize to the fractional form by applying its 
Maxwell formulation USj. This formulation consists of a spring (obeying 
Hooke's law) and a dashpot (obeying Newton's law of viscosity) in series; 
this arrangement shows a simple spatial separation of the solid (elastic) and 
the fluid (viscous) aspects and it is too speciflc to describe real viscoelastic 
materials. However, the hierarchical arrangement of a number (in general 
inflnite) springs and dashpots is already sufficient jl8j . Note, that in our ap- 
proach the spring defines a purely emotional or irrational investors' activity^ 
while the dashpot defines a purely rational one. 

* Usually the stress is denoted by a and strain by e. 

^ In psychology more often is used terminology 'affected driven activity' or 'authomatic 
activity'. 
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2.2. Conjecture 

There are several definitions of fractional differentiation and integration 
|19j . In what follows we are dealing strictly with the Liouville-Riemann (LR) 
fractional calculus. The fractional integration of arbitrary order a(> 0) of 
function f{t) is a straightforward generalisation of (jSJ, 

where qD^"' is the LR fractional integral operator of order a P^ . 

The fractional generalization of eq. ((T)) is performed by replacing expres- 
sions Tq-^ oDt^U{t) and rf ^ o by Tq"" q A~"f^(i) and rf ° o A""^(*) 
ones, respectively, where the fractional (in general) exponent a is a free but 
most important shape parameter. Hence, we obtain the fractional integral 
equation which is able to describe both independent paths (the rising and 
relaxation) of local temporary peaks of WIG: 

X(y) - X(0) = -rf " oDy'^Xiy) - A' sgn{C) Tq-" o^, "f/(y) 

+ J^[U{y)-Um, (11) 

where it was tacitly assumed that sgn{C') = —1, while the independent 
variable 

f tMAX - 1 for rising: t < tuAX 
\ t — tMAX for relaxation: t > t^jAX- 

As both paths of any peak are assumed as independent ones we consider all 
parameters present in ea. Hllf) as (in general) different for different paths. 

For relaxation the first term on the rhs of ea. ()ll|) describes feedback 
where the retarded value of index infiuences the present one; this value is 
sensitive here to the past one due to the algebraic, integral kernel. The 
second term on the rhs of ea. (|ll|l desribes explicitly a financial market 
retarded influence on the index (or the stock price); the third term gives 
the instantaneous influence. However, for the rising path the situation is 
more complicated and ea. (|ll() constitutes only a formal, convenient way to 
describe it. As we prove, the flrst and third terms constitute mainly the 
basis of a dynamical structure of the local (temporal) maximum of WIG. 

2. 3. Relaxation and rising fractional differential equations 

To make a step towards the interpretation, we deflne the fractional dif- 
ferentiation of order 7(> 0) 



oDVit) = ^ ( oDrVit)) , (13) 
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which is considered to be composed of a fractional integration of the order 

a n — 7 (—1 < 7 — n < 0) followed by an ordinary differentiation of 
order n. Now, by ordinary differentiation of the first order of ea- Hllj) and 
by applying definition H13() we obtain the linear inhomogeneous fractional 
differential equation 

"^^^y^ -- -rf " oDl-'^Xiy) - A' sgniC) r^^ ,Dl'^U{y) 



dy 



+ A' (14) 
dy 



which describes well both paths of the studied peaks. 



2.3.1. Free fractional relaxation: the reference case 

We found that the well developed local maxima of the index can be fitted 
(except for the left-hand side of the first one and up to their oscillations and 
fluctuations) by an intermediate part of the ML function. In our case we 
obtained several values of exponent a for WIG's maxima and almost all 
of them (except one) are smaller than 0.5; note that the left-hand side of 
the first maximum is well fitted by the usual exponential function (or by 
assuming a = 1 in the MF function). 

In other words, the relaxation of almost all WIG local maxima can be de- 
scribed by the fractional relaxation equation by setting in ea. ()14|) coefficient 
A' = 0. Such a simplified equation is, of course, a fractional generalization 
of the standard relaxation equation whose solution has indeed the form Q 
(where variable t is replaced by y and parameter r by corresponding ri one). 
This solution is considered here only as a reference case. 



2.4- Full Solution of the Fractional Initial Value Problem 
We solve the fractional initial value problem ()14|) by assuming that 

U{v) = ^ [exp(z (w - Aw) y) + exp(z (w + Atj) y)] , (15) 

(where parameters [/(O), w, Alij > 0) and by applying the Laplace transform 
of a fractional integral operator. Namely, the Laplace transformation of 
eq. (fTT|) yields 



where the Laplace transform of the LR fractional integral operator was ap- 
plied here. By introducing the Laplace transform of H15|) into ea. (|16() and by 
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Table 1 . First set of parameters describing the temporal maximum A 



Process 


a 


n [td] 


tuAX [td] 


Lo [td-'] 


Auj [td-'] 


Rising 


1 


58 ± 1 


320 ± 1 






Relaxation 


0.41 


176 ± 1 


316 ± 1 


0.18 


0.009 


Table 2. First set of parameters describing the temporal maximum a 


Process 


a 


n [td] 


tMAX [td] 


to [td~'] 


Alu [td-'] 


Rising 


0.31 


721.5 ± 1 


241 ± 1 


0.063 


0.005 


Relaxation 


0.25 


6830 ± 10 


245 ± 1 


0.26 


0.025 



using the inverse Laplace transformation in the time domain we can obtain 
the real part of the solution. However, to compare the prediction of our 
model with empirical data it is sufficient to use only the lowest order terms 
in the exact solution, i.e. it is suficient to use the following approximate 
solution 



X{0)+A' C/(0) sgn{C) (- 



-I 1 + 



^' C/(0) sgn{C) cos(u; y) cos(Aa; y), 



(17) 



since the parameters w. Aw, which additionaly multiply the integral terms, 
were found to be so small that the integral terms are negligible. 



2.5. Comparison with empirical data and discussion 

In Figs l2l5l we compared the empirical data defining WIG's local maxima 
(denoted as A, a, B, and C in Fig^ with the predictions given by formula 
H17|) . The monotonic curves (obtained by using only the first term) present 
free solutions while the oscillating curves (obtained by using the whole ex- 
pression) the full ones, i.e. the free solutions decorated with mono-frequency 
oscillations (rising and falling paths of peaks a and B, respectively) or wig- 
gles (right-hand paths of peaks A and a as well as left-hand path of peak B). 
The values of the key parameters which we obtained are shown in Tables 
[TEP . 

In FigElare presented: (i) the Mittag-Leffler function fitted, for example, 
to the left-hand path of the empirical maximum B (this is the free solution 
taken, in fact, from Fig0J), and corresponding (ii) KWW law (lower curve) 



Note that numbers shown without errors mean that errors are neghgibly smalL 
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Fig. 2. The temporal, local maximum A (cf. Fig^J: the evolution of the daily 
closing price of WIG in the semi-logarithmic scale (dots) from the beginning of the 
WSE activity to the end of the first maximum on March 30, 1995; the presented 
time series consists of data points for 545 trading days. The tangent solid line -was 
fitted to the left-hand slope of the empirical maximum. The curve fitted to the 
right-hand slope -was obtained by using ea. (|17|l . The values of the corresponding 
key parameters are presented in Tabled 



Table 3. First set of parameters describing the temporal maximum B 



Process 


a 


n [td] 


thlAX [td] 


CJ [td~'] 


Aw [td-^] 


Rising 


0.40 


164 ± 1 


401 ± 1 


0.035 


0.005 


Relaxation 


0.39 


366.5 ± 1 


370 ± 1 


0.055 





as well as (iii) the Nutting law (upper curve). This is a typical situation 
valid both for the left- and right-hand paths of all considered peaks. It is 
characteristic that none of the peaks reached the fully developed scaling 
region of the return per unit time ^'"^^^^ (cf. section [1.1(1 . In this sense 
the considered peaks have a precritical character. 

There are several other features common for all peaks which should be 



Table 4. First set of parameters describing the temporal maximum C 



Process 


a 


Ti \td\ 


tuAX [td] 


Rising 


0.40 


249.2 ± 1 


608 ± 1 
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Fig. 3. The local maximum a (cf. Fig^: the evolution of the daily closing value 
of WIG in the semi-logarithmic scale (dots) dated from the 770*'' trading day 
(the beginning of this local maximum) to the 1130*'* assumed as the end of this 
maximum range; the time series consists of data points for 360 trading days. The 
curves fitted to left-hand and right-hand paths of the maximum were obtained by 
using ea. H17|l : the corresponding key parameters are presented in Table |21 



noted: 

(i) Both paths of any peak can be considered as independent ones and 
the location of turning point (extremum) as a random one. 

(ii) The considered peaks are asymmetric since: 

• the exponent a, which characterizes the left-hand paths of any 
peak, is larger than the analogous one characterizing the right- 
hand one, 

• parameter ri (defining the time unit) of the left-hand path is 
smaller than the corresponding one for the right-hand path. 

(iii) The location of the extremal point tMAX given by expression ((T7)) is 
(in general) different for left- and right-hand paths of each peak. 

Moreover, a frequency modulation signal is necessary and/or the influence 
of the signal outside the maximum should be taken into account to describe 
the beginning of the left-hand paths of peaks a and B (cf. Figs 01 and [l]). 

Concluding, we suggest that our approach can rationally decrease the 
risk of investment on the stock market since it is able to warn the investors 
before the stock market reaches a critical region. 
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Fig. 4. The local maximum B (cf. Fig^: the evolution of the daily closing value of 
WIG in the semi-logarithmic scale (dots) dated from the 1380*'* trading day (the 
beginning of the local maximum) to the 2080*'' day assumed as the end of this 
maximum range; the presented time series consists of data points for 700 trading 
days. The curves fitted to both slopes of the empirical maximum were obtained by 
using ea.l|T7|l. The corresponding key parameters are presented in Table El 
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Fig. 5. The local maximum C (cf. Fig^l: the evolution of the daily closing value 
of WIG (dots) dated from the 2750*'' trading day (the be ginning of the local maxi- 
mum) to the 3382"'^ day assumed as the end of the maximum range; the presented 
time series consists of data points for 633 trading days. The curve fitted to the 
left-hand slope of the empirical maximum was obtained by using ea. ljlTII . The 
corresponding key parameters are presented in Tabled] 
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Fig. 6. The local WIG's maximum B (cf. Fig^): the analysis of the evolution of 
the daily closing value of the left-hand path presented in the semi-logarithmic scale 
(dots). The middle solid curve fitted to the left-hand slope of the empirical maxi- 
mum was taken from Fig^ while the lower and upper curves are the corresponding 
limits for the short time given by the streched exponential- or the KWW law 
and for asymptotic time the power- or the Nutting law is given by lO- 



We thank prof. Piotr Jaworski from the Institute of Mathematics at the 
Warsaw University for his helpful discussion. 
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